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Abstract 

One- and two-jet inclusive quantities in hadron collisions have already been 
calculated to next-to-leading order accuracy, using both the subtraction and the 
cone method. Since the one-loop corrections have recently been obtained for 
all five-parton amplitudes, three-jet inclusive quantities can also be predicted 
to next-to-leading order. The subtraction method presented in the literature is 
based on a systematic use of boost-invariant kinematical variables, and therefore 
its application to three-jet production is quite cumbersome. In this paper we re- 
analyze the subtraction method and point out the advantage of using angle and 
energy variables. This leads to simpler results and it has complete generality, 
extending its validity to n-jet production. The formalism is also applicable 
to n-jet production in e^e~ annihilation and in photon-hadron collisions. All 
the analytical results necessary to construct an efficient numerical program for 
next-to-leading order three-jet inclusive quantities in hadroproduction are given 
explicitly. As new analytical result, we also report the collinear limits of all the 
two-to-four processes. 
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1. Introduction 



The production of one or more jets is an important phenomenon in large transverse 
momentum reactions at hadron-hadron and photon-hadron colhders Jet produc- 
tion has large rates and rich final state structure, thus providing an excellent testing 
ground for the predictions of perturbative QCD and a formidable background to var- 
ious new physics signals. It is therefore important to have as significant theoretical 
predictions as possible for all the measured distributions in jet production. 

In the last decade the formalism of the next-to-leading order perturbative QCD has 
successfully been applied to explain quantitatively the data on one-jet and two-jet 
inclusive quantities obtained at SPS, TEVATRON and HERA At the next-to- 

leading order the theoretical uncertainties can be reduced below the experimental 
errors and impressive agreement has been found between the data and the theoretical 
predictions for all available inclusive one-jet or two-jet measurements. 

We recall, however, that there are two embarrassing unresolved discrepancies. 
First, the ramping run at TEVATRON could be used to test Feynman scaling; the 
data show larger scaling violation effects than predicted by the theory . Second, the 
tail of the transverse energy distribution of single-inclusive jet production, measured 
by the CDF collaboration 0, appears to be higher than the next-to-leading order 
QCD prediction (see, however, ref. |1T0|). 

The success (and also the possible difficulties) of the next-to-leading order de- 
scription of one- and two-jet production, and the large number of three-jet events 
collected by experimental collaborations motivate the extension of the quantita- 
tive next-to-leading order analysis also to the description of the three-jet data. The 
study of the ratio of three- to two-jet production rates is expected to give a clean 
measurement of as at hadron colliders. 

The theoretical analysis of inclusive three-jet production is rather complex. The 
five-parton amplitudes have to be calculated to next-to-leading order. As a result of 
significant technical developments based on the helicity method, supersymmetry and 
string theory, these calculations have been performed and presently the theoretical 
inputs at the amplitude level are fully available ||12| - p!5| . The leading order six-parton 
amplitudes, which have been obtained long time ago 0-|l8|, are also needed. The 
relation of the physical cross section with the partonic amplitudes, however, is rather 
involved, since the soft and collinear singularities appearing in the next-to-leading 
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order amplitude of the 2^3 parton scattering and in the leading-order 2 — 4 ampli- 
tudes are cancelled only when the cross section of inclusive, infrared-safe quantities is 
considered. In an efficient phenomenological numerical study, the soft and collinear 
singularities appearing in the loop corrections and in the real contributions have to 
be cancelled analytically. It has been pointed out in refs. |T9| - pl|] that the necessary 
analytical cancellation can always be achieved due to good universal properties of 
the soft and collinear limiting behavior of partonic amplitudes. In particular, two 



methods have been proposed: the so-called subtraction method [19,20| and the cone 



method ||2T|]. Both formalisms have successfully been applied to describe various 
distributions measured in inclusive one-jet and two-jet production. 

The aim of this paper is to present some further improvements of the formalism 
of the subtraction method necessary for the phenomenological study of three-jet in- 
clusive quantities. The subtraction method of ref. [|l9l is based on a systematic use 
of boost-invariant kinematical variables when the reference to the beam direction re- 
mains always explicit. As a result, the soft integrals become more complicated than 
their inherent structure would require and their treatment is somewhat cumbersome 
when it is applied to three or more jet production. 

We point out that the formalism of ref. [0] remains valid also when energy and 
angle variables are used. In this case the soft integrals become much simpler and 
the formalism is generally valid for inclusive production of any number of jets. The 
subtraction method with energy and angle variables has already been applied to the 
description of other production processes; for a formulation quite close to the one we 
use in this paper, see ref. P^ . 

The cancellation of leading singularities does not completely solve the numerical 
problems. The subtracted cross sections may still have integrable square-root singu- 
larities, thus resulting in an inefficient numerical evaluation. In ref. ||T^ this difficulty 
has been overcome by decomposing the cross section with partial fractioning into 
terms of single- singular factors. This method is not practical in the case of three-jet 
production since it leads to very long algebraic formulae. We shall show that such 
a decomposition is not actually necessary, since the measurement function defining 
infrared-safe inclusive jet quantities automatically splits the squared amplitude into 
terms where only one (or at most two) singular regions can contribute. 

In this paper we summarize the results obtained for the next-to-leading order five- 
parton amplitudes, we construct local subtraction terms for the singular regions of the 
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six-parton cross sections and we calculate the contributions of the soft and collinear 
regions analytically. With the latter result we have all the analytic formulae at our 
disposal to construct an efficient numerical program for the evaluation of three-jet 
production rates to next-to-leading order accuracy. 

Our paper is organized as follows. In section 2 the organization of the calculation is 
presented: our definitions and conventions are given, the MS collinear counterterms 
are explicitly worked out, the adopted jet-finding algorithm is described and the 
measurement functions are constructed. In section 3 we remind the known results for 
the virtual contribution. In section 4 we deal with the real contribution, and we show 
how to decompose it into single-singular factors. We define the soft subtraction using 
angle and energy variables and we evaluate the soft integrals analytically. We next 
turn to the definition of the subtraction terms for the initial and final state collinear 
singularities; all the needed integrals are analytically calculated. In section 5 we 
collect our results, and we show that when we sum the contributions of the real term, 
of the collinear counterterms and of the virtual term all the singularities cancel. The 
numerical implementation of our analytical results is also shortly discussed. Section 
6 contains our concluding remarks. The soft and collinear integrals are collected in 
appendix A. Finally, in appendix B we present a detailed description of the collinear 
limits of the n-parton cross sections. 



2. The jet cross section 
2.1. Introductory remarks 



Thanks to the factorization theorem [^, a generic differential cross section in hadronic 



collisions can be written in the following way 

da^H.H,)^^^^^^^^ = ^1 dx^dx2fi'''\x,)fi'''\x2)daMiKi,^2K2) , (2.1) 

ab 

where Hi and H2 are the incoming hadrons, Ki and K2 their momentum, and the 
sum runs over all the parton fiavours which give a non-trivial contribution. The 
quantities d&ab are the subtracted partonic cross sections, in which the singularities 
due to collinear emission of massless partons from the incoming partons have been 
cancelled by some suitable counterterms. 

We first express the subtracted cross sections in terms of the unsubtracted ones. 
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which can be directly calculated in perturbative QCD. To this end, we have to write 
the collinear counterterms. Due to universality, eq. applies also when the in- 

coming hadrons are formally substituted with partons. In this case, we are also able 
to evaluate the partonic densities, which at the next-to-leading order read 

/i'^ (x) = 6a,6{l -x)-^ [^Pa,{x, 0) - Ka,{x)^ + O (a^) , (2.2) 

where Pad{x, 0) are the Altarelli-Parisi kernels in four dimensions (since we will usually 
work in 4 — 2e dimensions, the in the argument of Pad stands for e = 0) and the 
functions Kad depend upon the subtraction scheme in which the calculation is carried 
out. For MS, Kad = 0. Writing the perturbative expansion of the unsubtracted and 
subtracted partonic cross sections at next-to-leading order as 

daa, = daiS + dal^^ , daa, = daf^ + , (2.3) 

where the superscript (1) denotes the leading (next-to-leading) order contribution, 
we have 

da^^^{h,h) = da^;^^{h,h) (2.4) 
da^^^{ki, k2) = da^Ji^ih, h) + J dx (^^Pdaix, 0) - Kda{x)^ da^^^{xki, fca) 

+ dx {^Pdb{x,Q) - Kdb{x)^ da^^Xki,xk2) . (2.5) 

The second and the third term in the RHS of eq. (|2.5| ) are the collinear counterterms 
we were looking for. Notice that in this equation the Born terms da^'^'' are evaluated 
in 4 — 2e dimensions. 

For three-jet production, the leading-order cross section can get contributions 
only from the two-to-three partonic subprocesses. We write this contribution in the 
following way 

X S3{{ki}3,5AJi}i,3)dMh,k2^ {ki}^,,). (2.6) 



Here we denoted with 

{<^/}l,3 = {-'^i, J2, J3} 



(2.7) 
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the set of the four-momenta of the jets. In the following, we will almost always 
suppress the indication of the {J/}i,3 dependence. We have indicated with ki and 
respectively the momentum and the flavour of the parton number i involved in the 
process; by definition, partons 1 and 2 are the incoming ones. In the sum J2{ai}3 5 
every a/, with 3 < / < 5, takes the values g, u, u, and so on. To avoid overcounting 
in physical predictions, we inserted the factor 1/3!. To shorten as much as possible 
the notation, we have collectively indicated the momenta as 

{fcj,,, = {ki\i<l<j}. (2.8) 

It will also turn useful to define 

{fcjij-] ^ {fc.l ^ < / < J, ly^njy^p,..}. (2.9) 

The same notation will be used for flavours. The quantity S3 is the so-called measure- 
ment function, which defines the infrared-safe jet observables in terms of the momenta 
of the (unobservable) partons; we will describe it in more details in the following. 
_^(3,o) jg ^YiQ two-to-three leading-order transition amplitude squared, summed over 
final state and averaged over initial state color and spin degrees of freedom, and 
multiplied by the flux factor 

2ki ■ ^2 uj{ai)uj{a2) 



= ^ ^n^^ E \^''-\^ - 3) % (2.10) 



color 

spin 



where uj{a) is the number of color and spin degrees of freedom for the flavour a. We 
remember that, in 4 — 2e dimensions, 

^(g) = 2iV,, uj{g) = 2{l-e)D^, (2.11) 

where Dj^ = — 1 is the dimension of the adjoint representation of the color group 
SU{Nc). Since the incoming partons can play a very special role, we will also write 
the functional dependence of Ai^^'^^ in the following way 

M^'''\au a2, {ajs.s; h, k^, {hh,,) . (2.12) 

The transition amplitude for processes in which only gluons and quarks are involved 
is usually evaluated in the unphysical configuration in which all the particles are 
outgoing. The amplitude A^^'''^^\2 ^ 3) of eq. (CT) 

can be obtained from the 

amplitude 

^(*-)^Q ^ 5) ^ ^(*--)(ai, as, {a,}3,5; -k,, -k,, {ki}^,,) (2.13) 
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simply by crossing (for details on crossing, see appendix B); notice that eq. ( p. 131) is 
crossing invariant. Finally, in eq. ( p. 61) we denoted with d(f)^ the full (i.e., with the 5 
that enforces the conservation of four-momentum) three-body phase space. 

Coming to the next-to-leading order contribution, both the two-to-three and two- 
to-four partonic subprocesses have to be considered. As customary in perturbative 
QCD, we denote as virtual the contribution of the former, and as real the contribution 
of the latter: 



(2.14) 



with 



{«J3,5 



X 53({A;J3,5; {Ji}i,z)dUki, h ^ {kih,,) , (2.15) 



{«i}3,6 



(2.16) 



where 1S4 is the measurement function, analogous to S3, for four partons in the final 
state. A^*^^'^^ is due to the loop contribution to the two-to-three subprocesses 



2ki ■ k2 Lu{ai)uj{a2) 

X E 



color L 

Spin 



^{tree)^2 3) [A^^""p\2 3))* + (A^''''^\2 ^ 3)) * ^('""^^ (2 3) 

(2.17) 



while Ai^^^ is defined in terms of the two-to-four transition amplitude 



2/ci ■ /c2 uj{ai)uj{a2) 



(2.18) 



spin 



Notice that in eq. ( |2.16| ) the dependence upon the momentum and flavour of the 
additional parton was inserted. Also, the measurement function, as well as the com- 



binatorial factor 1/4!, had to be modified with respect to eq. (p.6|) and eq. (|2.15|) . 
The full four-body phase space was denoted with dcp^. 
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We can now go back to eq. (|2.5|) , to write explicitly the collinear counterterms for 
the three-jet production. Using eq. (p.6|) we get 



X M'^^'°\d, 02, {a;}3,5; xki, k2, {ki}3^5)S3d(f)3{xki, k2 {^Js.s) , 

(2.19) 

X M^^'°\ai,d,{ai}s,5;ki,xk2,{ki}3^5)S3d(p3{ki,xk2 ^ {ki}^^^) . 

(2.20) 

By construction, these quantities, when added to the unsubtracted three-jet partonic 
cross section, must cancel the collinear singularities coming from initial state emission. 



2.2. The jet-finding algorithm 

We now turn to the problem of the definition of the measurement functions S3 and 1S4. 
As a preliminary, we need some prescription defining the way in which unobservable 
partons are eventually merged into physical jets. To define a jet in terms of partons, 
it is customary to distinguish two separate steps: the clustering algorithm, which 
decides whether a given set of partons is mergeable into a jet, and the merging 
procedure, which defines the jet momentum as a function of the parton momenta. 
The clustering algorithm we choose to use was introduced by Ellis and Soper . It 



is a kj, algorithm specifically designed for hadron-hadron collisions. It is formulated 
in terms of the transverse momenta kiT and of the lego plot distances Rij of the final 
state partons 

d^ = kl, (2.21) 

Rij = {V^ - Vj? + {V^ - Vjf , (2.22) 

d,, = rmn{kl,k]^)^, (2.23) 

where the constant D is the jet-resolution parameter which value is set at convenience; 
in practice, 0.4 < D < 1.0. The algorithm is defined by means of an iterative 
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procedure. One starts with an empty list of jets and a list of protojets, the latter 
being in the first step by definition identical to the partons. Then, the quantities di 
and dij are evaluated for all the protojets and the minimum among them is found. If 
this minimum is di, then protojet i is moved from the list of protojets to the list of jets. 
Otherwise, if the minimum is dij, then protojets i and j are merged into a protojet. 
The four-momentum of the protojet is defined by means of the merging procedure 
of ref. it is the sum of the four-momenta of the two constituent protojets. The 
jet-finding procedure is repeated as long as there are protojets around. When the list 
of jets is completed, that is, the list of protojets is empty, all jets with kj. below a 
certain threshold are dropped. We choose the threshold to be a given fraction (which 
we denote by /) of the greatest squared transverse momentum; this means that the 
threshold has to be recomputed for every event. 

Other jet-finding algorithms are obviously possible P,^. Nevertheless, prelimi- 
nary studies indicate that, contrary to e~^e~ annihilation, in hadropro duct ion the 
prescription described here is favoured by the data. 



2.3. The measurement function 

Using the jet-finding algorithm of the previous section, we can explicitly define the 
measurement functions we need to construct the jet cross section. 

When there are only three partons in the final state, no merging is possible, and 
the partons themselves will eventually result in physical jets. Therefore 

^3 = ^5{j,k,l)9{min{dj,dk,di) - f max{dj,dk,di)) 

a{J) 

X e{Rjk - Dy{R^i - Dy{Rki - D^) , (2.24) 

where we introduced the shorthand notation for the 5 over four-momenta 

5{j, k, I) = 6{k^ - J,)6{kk - J2)6{k - J3). (2.25) 

The indices j, k and / take the values 3, 4, 5 and are different from each other. cr(J) 
denotes the permutation over the jet four-momenta { Ji}i,3. 

When four partons are present in the final state, the situation is somewhat more 
involved. To get three jets starting from four partons, only the following possibilities 
may occur: 
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• no merging, but one parton is dropped being below the hard scale (this contri- 
bution will be denoted by S^^^); 

• one merging, occurring in the first step of the algorithm {Slp); 

• one merging, occurring in the second step of the algorithm 

• one merging, occurring in the third step of the algorithm {Sl^^). 
Consistently, we will then write 



1,3 

i<j 



where 

S^'^ = j:^ij^k,l)J^^\ (2.27) 

a{J) 

4^ = j:S{i+j,k,l)T^t'\ (2-28) 

a{J) 

= j:j:^iiJ^k + l)J^^\ (2.29) 

k<l 

S?^ = EE'^(^>J>^- + 0^^f , (2.30) 

k<l 

and the indices j, k, and / can take the values 3, 4, 5, 6 being different from each 
other and from i 

{j,k,l} = {3, 4, 5, 6}\0}. (2.31) 
To write explicitly the T functions we introduce the following shorthand definitions 

m{..,i, .., j +p, ..,mn, ..) = min{..,di, ..,dj+p, ..,dran, ••) , 
M{..,i, ..,j+p, ..,mn, ..) = max{..,di, ..,dj+p, ..,dmn, ■■) , (2.32) 
and, consistently, we denote everywhere di = i, dij = ij. We also write for short 

m([a]) = m(i, .., j/c, ..) \/i ^ a, \/jk^a, (2.33) 

where da can be either da or dajs, and the minimum on the RHS is evaluated over the 
list of all the di and djk with the exclusion of da- Finally, we denote 

mia,..,K.}([c]), (2.34) 
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which is analogous to eq. (|2.33| ); the {a, ..,b, ..} inserted means that the Ust of di and 
dij over which the minimum is evaluated is such that it does not contain the indices 
a and b. We then have 



kl 



^(i,2) 



(i,3) 



jkl 



= e(m([z]) - z)e(/M(j, k, I) - i)9{m{j, k, I) - fM{j, k, I)) 

X e{R^k - Dy{R^i - OyiRki - D^) , 

= e{m{\i3]) - i3)e{m{i + j, k, I) - /M(z + j, A;, /)) 

X , - Z}2)^^(i?,+,v - D^)e{Rki - D^) , 

= 9{m{\i]) - t)9{m{z,j, k + I) - fM{z,j, k + /)) 

X e{m{,y{[kl]) - kl)e{Rk+i,, - D^) , 

= 9{m{\i]) - ^)9{m{^,J, k + I) - fM{z,j, k + /)) 

X eim^.^ilj]) - j)9iD^ - Rki) . 



(2.35) 



(2.36) 



(2.37) 



(2.38) 



2.4. Infrared singular regions 

It is well known that the perturbatively calculated QCD cross sections, even after the 
ultraviolet renormalization, have a divergent behaviour, arising from the regions in 
which a parton (either virtual, that is, exchanged in a loop, or real, that is, emitted 
and contributing to the final state kinematics) is soft or coUinear to another parton. 
In this section, we will investigate the behaviour of the measurement function in this 
regions (which we will denote as infrared singular regions) in the case when there 
are four partons in the final state. The singular regions eventually occurring are as 
follows 

• parton i {i = 3, 4, 5, 6) is soft (A;° 0); 

• parton i is collinear to the incoming parton 1 or 2 (z || 1 or z || 2); 

• parton i is collinear to parton j {i || j). 

From eqs. ( |2.35| )-( p.38| ), it is quite easy to prove that the JF functions are non vanishing 
only in the following singular regions 

^^f ^ ^°-0, ^||1, z||2, (2.39) 
^if''^ ^ - 0, ^0, z II J , (2.40) 
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{i,2) 



(i,3) 



none . 



none . 



Eq. ( p.40[ ) suggests the following decomposition 



c(l) 



s\]'^d{dj-di) + s\]^e{di-dj)- 



(2.41) 
(2.42) 

(2.43) 



in this way, the first term in the RHS of eq. (|2.43|) does not get contributions when 



j is soft, and the second one when i is soft. This in turn implies that, after some 
algebra, we can cast eq. (|2.26|) in the following form 



where 



S 



(sing) 



sf^ + Y.s\^^e{d,-d 



(sing) _|_ ^(fin) 



(1)/ 



S 



{fin) 



Sf) + C(3) 



(2.44) 



(2.45) 



(2.46) 



and the [i] in the sum means that j can take the values 3, 4, 5, 6 with the exclusion 
of i. The quantity sl^^^"* does not get any contribution from the singular regions; on 
the other hand, ^f^*"^^ is different from when i is soft (or i || 1, i || 2), but it is equal 
to zero when any other parton is soft (or collinear to the incoming partons), thanks 
to the factor 9{dj — di). The region in which i || j contributes to s'f^^^'^ and to j5j'^*"^\ 
but the 6{dj — di) inserted in eq. ( p.45|) prevents any double counting. 

We can finally investigate the form of the limiting behaviour of the measurement 
function ^4 in the infrared singular regions. 

• i is soft. From eqs. ( |2.39| )- (|2^4^) and eq. (|2.44| ), we have 



lim ^4 



lim 



(2.47) 



The limit of the functions can be very easily evaluated from their definition; it 
simply amounts to make the formal substitution i + j —^j wherever i + j appears. 
We still need to do some combinatorial algebra to get the limit of the full 1S4; after 
writing explicitly the terms contributing to the sum and using 



lim 



(2.48) 
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we get 



where 



lim ^4 



(2.49) 



= T.^U^kJ)e{m{j,k,l)-fM{j,k,l)) 



(2.50) 



Notice that this quantity is identical to the 1S3 function of eq. ( |2.24| ), but for the fact 
that the indices j, k and / can also take the value 6 (and the value i is excluded). 
Eq. ( |2.49 ) has an obvious physical meaning: when a parton gets soft, the remaining 
partons act as physical jets. The measurement function has then to coincide with the 
one defined for three partons in the final state. We also point out that each term in 
the sum in the RHS of eq. ( |2.47[ ) has a well defined soft limit: 



lim 5f ) 



lim S, 



(1) 



<S3{{i])(^{D'^ — Rij)d{Rik — Rij)0{Rii — Rij)- 



(2.51) 
(2.52) 



1. From eqs. ( p.39|) -( p^ and eq. ( p.44| ), we have 



lim S, = lim 5f ) . 

It is quite easy to prove that, in this limit, the following equation holds 

1 



and therefore 



lim 6(m([i]) — i) 

hi II fci 



lim54 = 53(b;]). 



(2.53) 



(2.54) 



(2.55) 



The case in which i || 2 is completely analogous and gives an identical result. 
• i II J. From eqs. and eq. (|2^ ), we have 



lira ^4 

ki\\kj 



lim SP 

k^ 1 1 kj 



(2.56) 



where we have used the fact that S^j is symmetric in the exchange of i and j and 

9{di - dj) + 9{dj - di) = 1. (2.57) 
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In the limit at hand, dij — > and therefore 

hmO {m{[i j]) -ij) = 1, (2.58) 

while all the other 6 functions containing dij in the list over which the minimum is 
evaluated are zero. Therefore we have 

\im S^ = S,m), (2.59) 

where 

a(J) 

X e{Rpf, - Dy{Rpi - D^)e{Rki - D^). (2.60) 

and the indices p, k and / in the sum can now take also the value 7, having defined 
kj = ki + kj (notice that this formal manipulation allows to maintain the notation 
used in the previous case; the meaning of eq.( |2.59| ) is that the final state collinear 
limit of the 1S4 function is again a jet-defining function of three partons into three jets, 
where one of the partons has four-momentum equal to the sum of the four-momenta 
of the partons becoming collinear). 

We finally point out that in the numerical implementation of the algorithm we will 
also need to consider the case where one parton is soft and collinear to an incoming 
parton or to a final state parton. In this case, we find 

lim lim^f = S^m, (2.61) 
hm lim Sff = S^m). (2.62) 



Notice that in eqs. ( |2.61| ) and ( p.62| ) the order in which the limits are taken is irrele- 



vant. Also, eq. (|2.61|) holds true when i || 2. 



In the following, we will use the properties of the measurement function 1S4 to 
disentangle the structure of the singularities in the real contribution. In spite of this 
fact, our method is completely general. In fact, we will basically rely only upon 
eqs. ( p.44 )-( 2.46 ), that is, on the decomposition of the measurement function into 



terms which get contribution from two singular regions at the worst. It should be 
clear that such a decomposition can always be performed, being essentially due to 
the infrared safeness requirement on the measurement function. 
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2.5. Organization of the calculation 

Given the definitions of tlie measurement functions and of the colhnear counterterms, 
the main problem we must tackle now is to arrange for all the singularities appearing 
in the perturbatively calculated cross sections to cancel analytically. Notice that this 
cancellation is achieved only after that the real and the virtual contributions are 
summed together, as in eq. ( |2.14| ). Also, it must take place already at the partonic 



level, that is, we do not need to take into account the convolution with the partonic 
densities in the intermediate steps of the calculation (to have a physical picture of this 
fact, one can imagine to have very peaked partonic distribution functions, becoming 
a 5(1 — x) in some limit). A possible choice for the reference frame in which the 
calculation is carried out is that of the partonic center-of-mass one; this results in a 
simplification of the kinematics. We begin by considering the singular structure of 
the cross section; it is apparent from their definition that the measurement functions 
are boost invariant; therefore, proving the cancellation of the singularities in the 
partonic center-of-mass frame, ensures that the same cancellation holds true also in 
the hadronic center-of-mass frame. After the cancellation of singularities, we are left 
with a subtracted partonic cross section, analogous to ddab of eq. (|2.1| ), defined in the 
partonic center-of-mass frame instead of the hadronic center-of-mass frame. Clearly, 
the two are related by a longitudinal boost. The simplest way to perform this boost 
is by numerical methods (we stress that there is no loss of generality: in all practical 
cases, the integration in eq. ( |2.1| ) has to be carried out numerically). 



Since the real and virtual part of the cross section are both divergent, we have 
to write them in a form in which the structure of the divergencies is very clearly 
displayed. To this end, the most involved case is that of the real contribution. Taking 
into account eq. (|2.44 ), we split the real part of the cross section as follows 

da^'^ =y ida^''''^^ + da^f'""^) , (2.63) 

i 

where 

rf^lfS = E -M^'n{«^}i,6)5f"^#4, (2.64) 

{"i}3,6 



^^13 = E M^'\{a,U)St'''^d<P,. (2.65) 

{";}3,6 

c^clflji'i in eq. (|2.64|) is finite, and it is already suited for numerical computations. On 



the other hand, all the singular contributions are contained in da^a^^S, eq. ( |2.65| ). 
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The form of ^f*™^^ implies that dai^TJ!} has the following singularities: i soft, i \\ 1, 
i II 2 and i \\ j, Vj ^ i. From the previous discussion, we know that the singularities 
due to initial state coUinear emission are cancelled by the addition of the collinear 
counterterms, eqs. ( |2.19| ) and ( [^.201 ). It is then useful to define a subtracted real 



contribution in the following way 

d^M =da^''^ + + d(T('="*-) (2.66) 

a-^a^ a-^a2 a-^a2 a-^a^ V / 

that is 

da^"-^ = y (da^^'^'s) + (2.67) 

i 

where 

^^(^ms) = ^^{^^ng) ^ i^^(cnt+) ^ \^^{cnt-) _ .2.68) 

The key point here is that no parton but i can be soft. This amounts to a very clean 
disentangling of the soft regions, and allows to factor out immediately the pure soft 
singularities, as we will show later. 

In the following sections, we will first summarize the results for the virtual contri- 
bution, known from the literature. Next, we will discuss the separation of the soft and 
collinear singularities in the real contribution. We will then be able to analytically 
evaluate the structure of these singularities. Finally, we will collect the results in a 
form suitable for numerical evaluation. 

3. Virtual contribution 

Although the calculation of the loop corrections to the two-to-three partonic sub- 
processes is quite involved, and the result is rather complicated, the structure of the 
divergent terms is indeed very simple (see e.g. refs. |ll9| , |27| , |28| ) 



where 



\^ n=l ^ n=l / 



+^ E log^^^A<(;:„°)({aJ,5;{^.}i,5) 

ri,m— 1 V 

n^ra 

+-^LYH{«Ji,5; {fcji,5) • (3.2) 
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In these equations, fi is the renormahzation scale and Q is an arbitrary mass scale, 



introduced by Ellis and Sexton in ref. [28| to facilitate the writing of the result 



The quantities A^^^'^^ and A4^^'^^ were defined respectively in eqs. ( |2.171 ) and ( |2.1CI|) ; 
the Mmn^ usually denoted as color-linked Born squared amplitudes. They are 
symmetric in m, n and satisfy the identity 

Mt'^\{al}^,5■,{kl}l,5) = 16C{am)n'M^'''\{al},,,■{kl}^,,)■, (3.3) 

n = l 

the reason for the I/Stt^ normalization factor (that can be freely chosen, provided 
that eq. ( p.3|) is suitably modified) inserted in front of these terms in eq. ( p.2|) will 
become clear in the following. The factor 

{Any r(l + e)r(l-e)2 



rri-e) rfi-2e) 



C(e^) (3.4) 

naturally results from loop integration. In eq. ( p.2|) we also made use of the flavour 
dependent quantities C(a„) and 7(a„); for SU{Nc) color group, they are 

Cig) = C. = N^, 7(.) = i^^^^^^ , (3.5) 

6 

C{q) =C, = , 7(g) = ^C, , (3.6) 

where Tp = 1/2 and Nf is the number of quark flavours. All the non-divergent terms 
in eq. ( ^.21 ) were collected in m!"^s \ The explicit results can be found in ref. |]T2[ for 



5g process, ref. [|T3[ for 4qlg process, ref. |l^ , p!5[| for 3g2q process; the calculations 
were carried out using helicity amplitude methods in the dimensional reduction (DR) 
scheme. In the DR scheme, and AA^mn evaluated in four dimensions, 

while in the conventional dimensional regularization scheme (CDR) they are in 4 — 2e 



dimensions. As pointed out in ref. to convert the results obtained in the DR 
scheme |T2|-|T3| into the CDR scheme, one has to use A^^^''^^ and M.^^^^ in 4 — 2e 
dimensions instead of four dimensions; furthermore, the finite part is modified as 
follows ^ 

A^iYHcDR) = -MiY^DR) - E 7(an) , (3.7) 

n=l 

where the universal factors 7 are 

7to) = f, = (3.8) 



4. Real contribution 



4.1. Separation of singularities 

We now turn to the real contribution, defined in eq. ( p.l6| ) and further decomposed 
in eq. ( 2.63|) . Thanks to the fact that no divergencies are present in eq. ( 2.64 ), we will 



deal only with eq. ( 2.65| ). In the partonic center-of-mass frame, the incoming partons 
have momentum 

h = ^(1,0,1), (4.1) 



k2 = ^(1,0,-1), (4.2) 

where is the partonic center-of-mass energy and is the null vector in a (2 — 2e)- 
dimensional space. In this frame, we write the momentum of parton i as 



ki = U, \jl-yfeiT, yA , (4.3) 

where is a unit vector in the (2 — 2e)-dimensional transverse momentum space, 
— 1 < < 1 and < < 1. By construction, when — > the parton i gets soft, and 
when ?/j — > ±1 it gets collinear to the incoming partons. With this parametrization, 
in 4 — 2e dimensions the invariant measure over the variables of parton i is 

<(''''-k _ 1 (VS\ .1^2,1, ..2\- J, J. .n(2-20 



^^<''=(2^ri%=2(2^iVj i'^-y') di4y^-"\ (4.4) 

where dVlf ^^"^ is the angular measure in 2 — 2e dimensions. We also write the four- 
body phase space as (i</)4 = (i</)(i</)(i), where 

= (2.)--**-- + *2 - E n pIS^f ■ 

We now regulate the soft singularities by multiplying the invariant amplitude squared 
by Taking into account eq. (^^), eq. ( p.65| ) becomes 

= I ^ ^(^)({a,},6)5r^)rf0 



X 



2(27r; 



3-2e 



ii^'^' (l - yt) d^,dy,dQf'' ■ (4.6) 
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We can make use of the following identity 



-l-2e 



f-2e 
Scut 



(4.7) 



where ^cut is an arbitrary parameter satisfying the condition < ^cut < 1, and the 
distributions in the RHS are defined as follows 



< 



< 



/> 
/> 



6 



(4.8) 
(4.9) 



Notice that different E,cut can be chosen for different here we restricted to the simplest 
case. Substituting eq. ( |4.7| ) into eq. ( [4.6| ) we get 



where 



X 



1 



(4.10) 



{"i}3,6 
2-2e 



2(2vr; 



3-2e 



2e 



'log 6 



(2-2e) 



(4.11) 



2-2e 



2(27r; 



3-2e 



{«;}3,6 



(2-2£) 



1 - 2/J d^idyidn 



(4.12) 



Eq. ( [4.11|) does contain the singularities due to parton i becoming soft, while eq. (|4.12| ) 
is free of them, but does contain collinear singularities. We can also fully disentangle 
the collinear singularities of daa^al,i- Using eq. (|2.39|) and eq. (|2.40| ) we notice that 
a single term in the 5^**"'^^ function gets contributions from one singular collinear 
region at the worst. Therefore, using eq. ( |2.45D , we can naturally split dai^al^i into 
several terms: 



(4.13) 
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where 



{«i}3,6 



X 



2-2e 



1112 ''J 



2(2vr; 



1 



3-2e 



1 - yf^ ' d^idyidQ 



(2-2e) 



(4.14) 



X 



2-2e 



1 - 2/J d^idyidVt 



(2-2e) 



(4.15) 



2(27r)3-2^ \ 2 ; 

By construction, eq. ([4.14|) contains only initial state collinear singularities, while 



eq. (^4.151) contains only final state ones. 

In summary, the decompositions defined in eqs. ( p.63| ), ( [4.10| ) and ([4.13|) allow to 
write the real contribution to the partonic cross section in the following way 



a-^a2 / I a-]^a2,i a-]^a2,* a-]^a2,i / j a^a2,ij 



(4.16) 



where the quantities in the RHS of this equation have been presented in eqs. (p.64|) , 
( |4.11| ), ( |4.14| ), and ( |4.15| ) respectively. As far as the real contribution is concerned, 
eq. ([4.16|) will be regarded as our master equation. The key point is that every term 
in the RHS gets contributions from two singular regions at the worst. In the next 
three subsections, we will calculate the singular contributions explicitly. At the very 
end, we will show that, as expected, the sum of the real contribution, of the virtual 
contribution and of the collinear counterterms is finite. 



4.2. Soft singularities 



We begin by considering eq. ( [4.11| ). The 5{^i) allows to take the soft limit of all the 
quantities appearing in that equation. In particular, we have 



lim>[W({a,}i,6;{fcJi,6)=5, 



gai 



27r ^ 

n , m 



kn-km A^(3,0) /r^ 1 [i] . JT, iH 



6 ' 



(4.17) 
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where n and m in the sum run also over the values 1 and 2 (incoming partons). In 
eq. ( [4.171 ) we inserted /i^*^ as customary in 4 — 2e dimensions. Eq. ( |4.17| ) also explains 
the reason for the normalization of the A^mn'* terms in eq. (|3.2|) ; we choose to keep 
the form of the soft limit as simple as possible. Taking into account eq. ( |2.49| ) and 
the formal limit 

^lim rf0 = rf03 (^1, ^2 {^ijs.e) , 
and using I]{aj3g 5ga^ = J2t^ , eq. ( [4.11| ) becomes therefore 



(4.18) 



-2e 
cut 



)2£ 



s 



27r 2e 2(27r)3-2^ 



4' 

n<in i'^iJ3.6 



X 5{i.) {^] \ (l - y^)-' didmd^r-' ■ (4.19) 

The dependence of the RHS of eq. ([4.19|) upon the variables t/i and Qi is fully contained 
in the last line. The integral over dC,i, dyi and dVtf ^^"^ can therefore be performed 
explicitly; the result is reported in appendix A. We can cast eq. ( [4.19| ) in the following 
form 



X 



^ E '53([^])#3(fci,^2-{A:z}!;;6), (4.20) 



with 



f-2e cy2e 



2e 2(27r)3-2^ 
where J'mn is given in eq. (|A.8|) . Explicitly 



o Urr. 



(4.21) 



'j-{div) 
ran 

jireg) 



1 (47r)- //.^ 

87r2r(i-e) l^g2 

1 

8^ 2 



^ - - (log 



log 



1 /■2 Q ;^2 C t . t 
log^ + log ^ log 1^ - Li, 



1 , 2 2fcri ■ km 

H — log 

2 E E 



log (4 



2,EjiEm 



E E 



log 



AEnEm 

kn ' km \ 

lEfiEm / 



h ■ h 

f^n f^m 

2E E 
^-'-'n-'-'m 



(4.22) 



2 log^ 2 



. (4.23) 



-21- 



Here En is the energy of the parton n in the partonic center-of-mass frame and Q 
is the mass scale introduced in eq. ( p.2| ) . We can express En in term of invariant 
quantities through the equation 



It is apparent that the dependence upon the index i in eq. ( [4.20| ) is immaterial. We 
can therefore relabel the partons at will and get 

5 

27r 2 



d^t. = E ^^l:L. = ZIT. ifi^ + ^i^O E da^lXi^iU) (4.25) 



where we have set 



d<7^1 ({aji.s; {ki}i,,) = ^-M^^f ({az}i,5; {k}i,5) dcl>, (h, k, ^ {kh,,) (4.26) 

(in the following, as in eq. ( [4.25| ), we will not indicate explicitly the dependence 
upon the momenta of this quantity). Exploiting eq. ( p. 3D and using the fact that 
E^ = E2 = y/S/2, eq. (g) becomes 



(,) ^ as {4nY 
"^"^ 27rr(l -e) lg2 1 ^ 



Ms. 



1 5 2 A . , 2Ei 



as (47r) 



271 Ti 



2E^K) + -E^(%)iog— ^ 

n=l ^ j=3 ^CMtV O 

c^cT(°H{aJi,5) 
47r)^ //i^Y 1 5 2A;„ ■ A;^ 1 



-(C(ai) + C(a2)) log ec«t 
e 



+ 1^ ^ E ^i^^ E da^liia,}^,,) . (4.27) 

n7=m 



where, consistently with eq. 

rfa(°)({««}i,5; {A;Ji,5) = {{aih,,; {k},,,) d^, (h, h ^ {hh,,) ; (4.28) 



this quantity is just the Born cross section of eq. ( p.6|) without the sum over the 
flavours of the final state partons. We point out that the second and the third 
singular term in eq. ( [4.27 ) are cancelled by similar terms appearing in the final and 



initial state contributions. The remaining singular terms are cancelled by the virtual 
contribution. 
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4.3. Initial state singularities 



We now turn to eq. ( [4.14| ). The divergencies due to parton i becoming coUinear to one 
of the incoming partons can simultaneously be regulated by multiplying the invariant 
amplitude squared by the factor (1 — yf). Eq. ([4.14|) becomes 



1 
4! 



X 



2e 



2-2e 



{"i}3,6 



2(27r)3-2^ \^ 2 
We use the identity 

-1- (25,)- 



{1-yn 



2e 



6{l-yi) + 6{l + yi) 



where 



r{y. 



1 



+ 



(4.29) 



(4.30) 



(4.31) 



Here 5/ is an arbitrary parameter satisfying the condition < 5j < 2, and we defined 
/ 1 ' 



< 



< 



/> 
/> 



^,^/(.,.) -/(!)%. -1 + 5.)^ (4.32) 
-1 1 + 2/i 



Very much like the ^cut parameter, it would be possible to choose a different 5j for 
every we prefer however to avoid the proliferation of such parameters. Eq. (|4.29| ) 
can be splitted into three terms, the first two which contain the coUinear singularities 
(and in which the invariant amplitude is substituted with its coUinear limit), and the 
third one which is finite in the limits yi ±1. Explicitly 



a-^^aQ^^i a^a2,'i a^a^,^ a-^a-^^t ' 



(4.34) 



where 



2-2e 



2(2vr; 



3-2e 



X 



(4.35) 
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and 



,6 



2e 



' log 6 



2-2e 



3-2e 



(4.36) 



Thanks to the subtraction prescriptions of V{yi), which regulate initial state diver- 
gencies, and to eq. (|2.39|) , which guarantees that S^^^ does not get contributions from 
final state collinear divergencies, eq. ( [4.35|) is finite, and therefore we can set e = 0. 
We get 

ii/i\r/ 1 \ / 1 \ 1 1 /,/c\^ 



a-i an ,1 



4! 2 ve. 



1/i 



1 + 



2(27r) 



(4.37) 



Since all the divergencies have been properly regulated, the RHS of eq. ( [4.37| ) can be 
numerically integrated. 

We still have to deal with the divergent part of eq. ([4.34|) , that is to say, with 
the quantities defined in eq. ( |4.36| ). The 5(1 =F yi) in eq. ( [4.36| ) allows to take the 
appropriate collinear limit of all the quantities appearing in that equation. We only 
consider daa^a^}, that is the case in which i || 1; the treatment of i || 2 is completely 
analogous. The collinear limit of the invariant amplitude squared can be written as 
(see appendix B) 

hm M^'\{ai},,e; {kU) = ({a,}i,6; {k},,,) 

+ ^)-^^''°^ «^)' zh, {h]%) , (4.38) 

where Pabi^-, e) is the Altarelli-Parisi kernel for z < 1 in 4 — 2e dimensions and S'(c, d) 
is the flavour of the parton which can split into two partons of flavour c and d\ for 
example, S{g,g) = g and S{q,g) = q; by definition, if the splitting into the flavours 
c and d is not possible, Ps(cd)a zero. The z parameter is such that, in the collinear 



limit, ki = (1 — z)ki. Taking into account eq. 



we have 



^i = l- z if yi = l 



(4.39) 
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Using eq. ( |4.38|) , we can then write 



(4.40) 



In this equation, the dependence of A upon the azimuthal variables of parton i is 
trivial; therefore, it can be shown that 



dfif ~^'^A = 0. 



(4.41) 



Furthermore, the second term in eq. ( |4.4CI| ) does not depend upon the azimuthal 
variables, and after the dfij-^"^^^ integration in eq. ( |4.36| ) we get a factor 



(2-2e) 



2tt 



1-e 



r(i-e) 

Putting all the numerical factors together, taking into account the formal limit 



(4.42) 



lira 

ki\\ki 



using eq. ( p.55| ) and 



-fE + login 



(4.43) 



(4.44) 



we get 



E 

{«l}3,6 



, Si,- 



2e 



'log 6 



X 1a<(3.o) ^Sia,, a,), a2, {ai}%; (1 - i,)h, h, 



(4.45) 



When the sum over Oj is performed in eq. ( |4.45|) , the quantity S'(ai, a,) takes all the 
possible flavour values such that Ps{a^ a,)ai different from zero. Therefore, eq. ( |4.45| ) 



-25- 



can be rewritten as 



a 



ai ao ,t 



S5/ 



2e 



'log 6' 



e.P<^(l-e.,e) 



X 



i E -M(=^'°)(rf,a2,{a,}!;;e;(l-e.)A:i,A;2,{A;jgl 



(4.46) 



X 53([^])#3 ((1 - ^1, h {kiYi],) di. 

We can now show that the divergent part of eq. ( |4.46|) is cancelled by the coUinear 
counterterm of eq. ( |2.19| ) up to a factor 1/4 (see eq. ( |2.68| ) to understand the 1/4 
factor). To this purpose, we rewrite eq. ( p.l9| ) in terms of the variable ^j. Using 
eq. ([4.39|) , it is a simple matter of algebra to prove the following identity 

1 



5(l-y,,)(^) =5(l-y,) 



+ 5(e.)ioge 



cut 



(4.47) 



where the 5{1 — Ui) formally states the fact that the above equation is strictly valid 
only in the coUinear limit. The Altarelli-Parisi kernel can be put in the form 

(l-z)P<(;^,0) 



+ 7(0)5^6(^(1 - z) 



(4.48) 
(4.49) 



2C{a)5ab5{l - z) 

where the quantities C{a) and 7(a) were given in eq. ( ^.51 ) for gluons and eq. 
for quarks. Using eq. ( |4.47| ), ( |4.48| ) and ( [4.49| ), and changing the integration variable 
2; — s> 1 — ^j, it is straightforward to see that eq. (|2.19|) can be cast in the form 

d \ 

^^Pda,{^ - 0) + 5 ,aAid {l {d) + 2C{d) logC 



1 

+ - 

e 



X S;{\i])dh ((1 - ii)ki. h ^ {k}f^ di, , 



(4.50) 



since the dependence upon the index i is immaterial. Eq. ( [4. 501) is of the same form 
of eq. ( |4.46| ). It is now apparent that the quantity 



(4.51) 
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does not contain any purely collinear pole (i.e., proportional to the Altarelli-Parisi 
kernel for z < 1). Expanding in a Taylor series eq. ( |4.51| ), and using 

Pati^, = P<{z, 0) + eP^<(z, 0) , (4.52) 

we have 

a. 1 



da^X^ = g-(7(ai) + 2C(aOlogec 



+ ^EU.p.t, (1-6,0) 



'log 6 



-eii'i<(l-6,0)(^-j -AV(l-e.)| 



(4.53) 



xSsiltm, ((1 - ^,)h, h ^ {kK]e) • 

The pole part of this equation, which originates from the 6 term in the flavour-diagonal 
Altarelli-Parisi kernels, cancels a corresponding term in the soft-virtual contribution. 
The remaining part is finite, and can be numerically evaluated. The dependence upon 
the index i in the RHS of eq. ([4.53|) is trivial, and the sum over i is very easily carried 
out. After some algebra, we get 

m,+) 



i 



--] i^l 

27rVr(i-e) VQV e °^gv 

x(7(ai) + 2C(ai)logUt) E rf^^°H{««}i,5; {fc/}i,5 



a 



+ ^EU^.l.(i-e,o) 



1 



S5, 



log^ + 2 



'logC' 



-ePl<(i-60)(^-J -ir,„^(i-o| 

X ^ rfaW(^^,«2,{az}3,5;(l-0^i,^2,{A;z}3,5)t^e, (4.54) 

{ails, 5 
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where do"*^'^^ was defined in eq. ( [4.28| ). In the very same way, we have 

27T\T{l-e) [qy e ^ ) 



+ |E{eP.<.(i-e,o) 



X E f^a(°)(ai,t^,{««}3,5;A;i,(l-0^2,{A;,}3,5)f^e- (4.55) 

{«;}3,5 



4.4. Final state singularities 



We finally turn to the problem of regulating the divergencies appearing in eq. (^4.15|) . 
To this purpose, we rewrite eq. ([4.3|) in the following way 

h = (l, A;,) , k=vR, P= (O, l) . (4.56) 

Here p is a 3 — 2e dimensional vector and i? a 3 — 2e dimensional rotation matrix. 
With e = and eq. ( [4.3|) , the explicit form of R can be worked out, but we will not 
need it in the following. We then parametrize the momentum of parton j as 

kj = —ij (l, kj) , kj = pjR , pj = [^Jl-y]ejr, Vj) ■ (4.57) 
From this definition we get 

h-k,= (^y!pj e.O(l-%). (4.58) 

Therefore, the collinear limit i || j is obtained with the parametrization of eq. ( |4.57| ) by 
letting Uj 1. It follows that we can regulate final state divergencies by multiplying 
the invariant amplitude squared by the factor {1 — yj). We also write 

d(t) = d(f>d^{j) (4.59) 

2-2e 

^]-'^ (l - y') d^.dyjdQf-''^ , (4.60) 



S 



2(27r 
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where, from eq. (^]5 



/ 6 \ [ij] rl^-2eh. 



Eq. ( [4.15| ) becomes 



1 

4! 



2e 



{«i}3,6 



(4.61) 



X 



X 



_l 

2(27r)3-2^ 1^ 2 



2-2e\ 2 



We can therefore use the identity 



(4.62) 



(1-%) 



-l-e 



-6{l-yj) + 



+ 0{e) 



(4.63) 



where < (5o < 2, to spht (icri""*^^^ into two terms 

da^'^*\ . = (icT(°"*'+) + da^""*'^) , (4.64) 

where dai°al'tj^ is proportional to 6{1 — yj) and dai°^l\{] singularities (and 

therefore we can set e = 0) and can be numerically integrated. Explicitly 



1102 ,ij 



X 



and 



4! 

1 /v^ 

2(27r)3 1^ 2 



1112, »J 



5o 
2\ 2 



{«i}3,6 



d(j)dC,id^jdyidyjd(pid(pj , 
{26 oY 



-5(1 — yj)dyj VA dfj, . 



where 



= -r 
4! 



'log 6 



C - 



{ci!}3,6 



(4.65) 

(4.66) 

(4.67) 
(4.68) 
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2(27r)3"2. I 2 



d^]-^^ (l - y^) ' d^d^.dydnf^dn 



(2-2e)^(^(2-2e) 
j 



(4.69) 
(4.70) 



We perform the following change of variables 

6 = (1 - z)^7, = z^7, diidij = ^rd^jdz . 
With this definition, we have 

\/~S / \ 
5(1 - y,) + k,) = 5(1 - y,)kr = 5(1 - %)^6 (l, ^l-yHr,y^) , (4.71) 

which defines the momentum k^; notice that the definition of ^7 is therefore consistent 
with eq. ( ^.3| ); in the collinear limit, this parameter is proportional to the energy of 
the parton which eventually splits into the two collinear partons i and j. With the 
definition of z in eq. ( |4.70|) we have, analogously to eq. (|4.38| ), 



+ 



k-i ' kj 



-A ({az}i,6; {ki} 



1,6) 



-^a%(^,e)-M^''°H{«jS:^;{A:.}S:^), (4.72) 



with 



07 = S{ai, aj). 
Using eqs. (gT^), ( P3E| ) and ( ^3^ we get 



(4.73) 



lim A = Anasfi 



2e 



l-z 



9 z 



since from eq. ( 4.7CI|) immediately follows 

e{d, -d,) = 9(^z-^y 

We can now exploit the following identities 

^7 



^7 



(4.74) 
(4.75) 

(4.76) 
(4.77) 
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where 



log(l - z] 



1 - z 



+ log 6 



l-z 



(4.78) 

^ (logger- log' Ui) 5(1 -^). 

(4.79) 



Finally, with the change of variables in eq. ( [4.7CI| ), the measure becomes 

2-2e\ 2 

/ I /^/,S'^ 



2(2vr; 



3-2e 



(4.80) 

After some algebra, and using the analogous of eqs. ( [4.41| ) and ( |4.42| ) when integrating 
over dQ^^ '^^\ we get 



e r(l -e) I 2 ; 4! 



p(°)(z) -26^(1) (; 



X rf.(l-.).-'^ E ^a%(^,e)-M^^'°H{«'}i^;{^.}r^)'53([^j]) 

{li}3,6 



1 



2-2e\ 2 



2(27r: 



3-2e 



^er^'^c^^T {l-y^y'dy,dQ. 



(2-2e) 



(4.81) 



Using eq. ( |4.71|) , we can see that the following relation holds 

2-2e 



dh {h,h^{ki}t'}) 



1 (Vs' 

2(27r)3-2^ 1^ 2 ^ 



(2-2e) 



(4.82) 



Therefore we have 



-'"ajar '^^-^ajaj 



where 



(0) 



X, 



(1) 



^A^^^^'") ({a,}f^ ; {A:,}f,^) 53([u])#3 (A^i, fc, - , (4.83) 

= f^dzz-'^{l-z)P<{z,e)e[z-]^V^'\z), (4.84) 
= /^'rf^^-'Ml-^)i^at(^,e)^(^-^)^^^'H^)- (4.85) 
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Eq. (|4.83|) shows that the singular term due to final state coUinear emission factorizes 
into the product of two quantities; a cross section, which has a two-to-three partonic 
kinematics, and a term dependent upon the Altarelli-Parisi kernels, which completely 
describes the collinear splitting. We can sum over i,j the quantities dai°^l'P . ^ 
simple matter of combinatorial calculus to see that this sum can be rewritten in the 
following way, by simply relabeling the partons 



i=3 



(4.86) 



where 



a 



« 3 



S5o' 



(1) 



rfaW({aJi,5). (4.87) 



^ {a-l}3,5 d 

The expression for Y.d^d'aj ~ 2eX^^|] is given in appendix A. By substituting it in 
eq. ([4.87|) we get 



as (47r)^ //i^ 
27rr(l-e) 



+ 



X ■ 



E 



7(«j) - 2C(aj) log 



2E, 



5 

xE 

i=3 



7'(aj) - log^ (7(0^) - 2C{aj) log 



2K- 



2g 



J2 



+2C(a,) l^log^ - log^ J - 27(a,) log 
The explicit expression for 7'(aj) is given in eqs. (|A. 12|) and (|A.13|). 



2Ei 



5. Results 

In the previous sections we have shown that the subtracted next-to-leading order 
contribution, defined in eq. ( |2.5| ), can be written as the sum of a term with a two-to- 
four partonic kinematics, and a term with a two-to-three partonic kinematics 

a-]^a2 a.]^a2 '^l'^2 
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The two-to-four part is 

^^(i,4p) = y ( da^f'^'X + ticr^'"'^) + y , (5.2) 

where the quantities appearing in the RHS of this equation were given in eqs. 
( |4.37|) and (f4.65|) respectively. The two-to-three part is 



da(^'^P) = rfa(^) + da^'^ + + rfa^*"'") + c/cr(°"*'+) (5.3) 

a-^a^^ a-^a^^ a-^a^^ \ / 

where the quantities appearing in the RHS of this equation were given in eqs. ( |2.15| ) 
(with the matrix element of eq. (|]1]))' (EM), (^5^ ), ( F55|) and (E88|). Every term 



in the RHS of eq. ( p.2|) is finite, and can be numerically integrated. On the other 
hand, the quantities in the RHS of eq. ( |5.3| ) are divergent; nevertheless, from their 
explicit expression previously reported, it is apparent that the divergencies cancel 
in the sum, and da^}^^'^"* is finite. The finite part of eqs. ( 4.54 ) and ( 4.55|) , due to 



initial state collinear contribution, is slightly more complicated than the usual two- 
to-three kinematics contribution, since it has an additional folding in the variable ^. 
Therefore, we further split eq. ( ^.3|) into a sum of two terms, one of which does not 
contain any folding. The reason for doing so is that in a numerical computation the 
two terms have to be differently treated, in spite of the fact of having the same final 
state kinematics. We have 

cia(^'3p)^^^(i,3H^^^(i,3pO. (5.4) 

a^a2 a-^a2 a-^a^ V / 

The part without ^ folding is 

^^''^ = S ^ Q({aJi,5)rfaW({a,}i,5) 

{a; Is, 5 

+ E -M(,YH{«Ji,5)53#3, (5.5) 

Wt3,5 
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where 

2({ai}i,5) 



E 

j=3 



2Q2 



7(aj) — 2C{aj) log 



+2C(a,) log 



2 



5 



log ^cut - 27(aj) log 



2E,' 



Vs 



- hg^i^{ai) + 2C{ai)\ogUt + l{a2) + 2C{a2) hguX (5.6) 



and da^^\ da^^j^ were defined in eqs. ( |4.2(j| ) and ( ^.261 ) respectively. The part with the 
folding can directly be read from eqs. ( 4.54 ) and ( 4.55| ), and it is 



'log^' 



-ep^<(i-e,o)(^^^ -ir,,,(i-o} 

X J2 d(J^°\d,a2,{ai}-i^5;{l-^)ki,k2,{ki}-i^5)d^ 



+ ^EU^.a.(l-e,0) 



'1\ , 5(5, 



'log^' 



-ep^<(i-e,o)(^-j -K,,,(i-o| 

X da^^\(^i^d,{ai}-i^5;ki,{l-^)k2,{ki}-i^5)d^. 



(5.7) 



We have derived our results in the partonic center-of-mass frame. This slightly 
constrains the general validity of the formalism, since in the numerical computations 
Xi and X2, the Bjorken x of the incoming partons, have to be chosen as independent 
integration variables. Notice that they define the boost from the hadronic center-of- 
mass frame to the partonic one. We stress that this constraint can be very easily 
relaxed, without affecting the correctness of the derivation. Nevertheless, as far as 
the numerical computations are concerned, it appears to be advantageous to work 
in the partonic center-of-mass frame. We remind that in previous applications of 
the subtraction method to jet physics Xi and X2 have been expressed in terms of 
transverse momenta and rapidities. With our finite next-to-leading order partonic 
cross sections, we can get the physical cross section using eq. (|2.1|). 



-34- 



A final remark on the numerical computation: the parameters S,cuti and 5o 
define the soft and collinear subtractions. Thanks to this property, we can freely 
redefine the terms in the RHS of eq. ( |5.1| ) without affecting the sum. This gives us 
the possibility of a significant numerical check on the correctness of the calculation; 
the physical results have to be independent from the value of these parameters. Also, 
notice that a clever choice of the parameters results in saving computing time. If for 
example eq. ( |4.8|) is considered, it is apparent that in most cases the soft counterterm 
(/(O)) is not evaluated if ^cut is small. On the other hand, ^cut should not have a value 
too close to 0, because in this way the quality of the convergence is rather poor. 



6. Conclusions 

We have studied the production of three-jet inclusive quantities at the next-to- 
leading order in QCD, using the subtraction method. Comparing to previous treat- 



ments |T^, in which the one- or two-jet production is dealt with this method, we 
used a somewhat different approach in at least two important aspects. First of all, we 
used angle and energy variables, instead of transverse momenta and rapidities; this 
results in a simplification when integrating the eikonal factors, which appear when- 
ever one of the final state partons gets soft. As a second feature, we fully exploited 
the measurement function, which defines infrared-safe observables, to disentangle the 
soft and collinear singularities in the real contribution. In fact, the measurement 
function can always be cast as a sum of terms, each of them getting contributions 
only from the infrared singular regions associated with a given parton. At the very 
end, we have therefore splitted the cross section into single-singular contributions, 
without the drawback of decomposing the transition amplitude squared with partial 
fractioning. 

The analytical results obtained were organized in a form suited for numerical 
computations. We splitted the real plus virtual contribution in a sum of terms with 
different final state partonic kinematics. The possibility has been left to redefine 
them up to a finite piece, without affecting the sum, which is the only quantity being 
physically meaningful. 

Although we treated explicitly only three-jet production in hadron-hadron colli- 
sions, our formalism can be easily extended to n-jet production. Furthermore, we 
emphasize that it is also applicable to n-jet production in e+e~ annihilation and in 
photon-hadron collisions. 
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APPENDIX A: Soft and collinear integrals 

In this appendix, we collect the results for the integrals of the eikonal factors, 
which are substituted into eq. ( [4.19| ), and for the integrals of the Altarelli-Parisi 



kernels, with the prescription defined in eq. ( [4.84| ) and ([4.85|) . 



Soft integrals: given eq. ( ^4.19 ), we will consider the quantities 



iJ^nm — J EnmdoJi , (^•-'-) 

where 

du;, = {l - y^Y' dydnf\ (A.3) 
Using eq. ( [4.3|) we can rewrite 

l-COS^n^ 
(1 - COSt^ni) (1 - COSt'mi) 

where Oa/B is the angle between the directions of the three-momenta and kj3. The 
measure dui is the angular measure in 3 — 2e dimensions for the momentum ki] since 
also the eikonal factor is rotationally invariant, we can redefine at will the angular 
variables of the momentum ki. If we choose 

/cj = (.., sin (yj sin 6^, cosy? sin 6^, cos 6^) , (A. 5) 

then we have 

duji = 2^-'^'7T-' "^/^ ~ ^\ sin-^^ if sin-^' edcosOdif, (A.6) 

r(l-2e) ^ ^' ^ ' 

where an integration has been carried out over the variables upon which the eikonal 
factor can not depend. Finally, we decompose eq. ( [A.4| ) as follows 



_ 1 1 / 1 



1 — C0S6'„,- ' 2 V "™ 1 — COS 6'. 



+ m). (A.7) 

1 - cos 6rr,.i J 



The first term in the RHS of this equation will give, after integration, a collinear 
pole. The second term is regulated by means of the subtractions, and will contribute 
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a finite term to the integral. The final result is 

rri-e) 



I vr 



rfl-2e) 



X 



2^^7rr(l - 2e) 

(r(i-6)f 

-2e7r ( -Lis 



-^ + 41og2 + 2e(y-21og22j j +27rlog 



1 — cos Qr, 



1 — COS Qr. 



-\oi (2(l-cosft„ 



— -21og^2 



+2e7rlog 2(l + cose„,„) log 



1 — cos Qrirr 



(A.8) 



It is just a matter of trivial algebra to get from this equation the expressions of 
eqs. O) and (lOSi) . 

• Collinear integrals: the calculation of the integrals defined in eqs. ([4.84|) and ([4.85|) 
is a lengthy but trivial operation. Using eqs. ( [4.78|) and ( [4.79| ), and taking into ac- 
count that the relabeling of the partons amounts to the substitution ,^7 C,j in the 



expressions of I^J. and 1^^'. , and that by construction 



r(l) 



2^ 



we get 
where 



d 



X 



(0) 



dai 



2eX 



(1) 



dai 



(A.9) 



(A.IO) 



2E- 

2C{aj) log ■ ^ 



^cutVS 



7(«i) 



7'(a.) - 27(«.) log ^ + 2C(a,) [log' -jl - log^ i 



cut 



(A.ll) 



with 



l\9) 



67 
13 



27r^ 



23 



27r2 



(A.12) 
(A.13) 
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APPENDIX B: Collinear limits 



In section 4 we have stated that the A term of eqs. ( ^.381 ) and ( [4.72| ) does not 



contribute to the result, since its average over the azimuthal angle of the collinearly 
emitted parton is zero. This fact can be proved in 4 — 2e dimensions studying the 
collinear limit in terms of the Sudakov variables (see e.g. ref. [|T5). In spite of the fact 
that the contribution of A is zero, it is of some importance to know its explicit form, 
since it appears in local subtraction terms where it improves the numerical treatment 
of equations like, for example, eq. ( |4.37| ). The goal of this appendix is to express A 
as a function of the helicity amplitudes known from the literature. We will work in 4 
dimensions, since we will eventually use A only in numerical computations. 

In the following, we will consider the production process of n particles 

1 + 2 ^ X + z + J (B.l) 

in the limit in which the two massless particles i and j become collinear to each other 
or collinear to the incoming massless particles 1 or 2. Notice that the properties of 
the n — 2 particles denoted collectively by X are unspecified and of no importance in 
what follows. This implies that the results of this section are completely general and 
their validity is not restricted to jet physics. 

We begin by considering the final state emission, that is the case in which i and 
j are collinear. Therefore, we study the structure of the process 

P ^ ij. (B.2) 

We have to consider only three cases: 

a) g ^ gg, (B.3) 

b) ^ ^ qq, (B.4) 
c) q ^ qg, (B.5) 

where in case c) q can be either a quark or an antiquark. We begin with case a). 
Giving color indices b and c to the splitted gluons i and j, and defining the fraction 
of momentum z through the equations 

ki = zkp , kj = (1 — z)kp , (B.6) 
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{h, 


ha, 


hb) 


9 ^ 99 


g ^ qq 


q ^ qg 


(+ 


+ , 


+) 


[ab] 







(+ 


+ , 




-z^ < ab > 


^1/2 (1 _ ^)3/2 < > 


-2;^/^ < ab > 


(+ 




+) 


-(1 - z)"^ < ab> 


-^3/2(l_^)l/2 





(+ 

















+ , 


+) 













+ , 




(1 - zf [ab] 




^3/2 (1 _ ^)l/2 









+) 


[ab] 


-^1/2 (1 _ ^)3/2 


z3/2 [ab] 








- < ab> 





-z^l'^ < ab > 



Table 1: Splitting process P{h) a{ha)b{hh) for all the possible choices 
of the partons a and b and the helicities h, ha and hb- The splitting func- 
tions S^p"^'' are obtained by dividing the entries of the table by < ab > 
[ab]^y z(l — z). By construction, the parton a has the fraction z of the mo- 
mentum of the parton P. The splitting functions Sgq can be obtained from Sqq 
with the formal substitution z ^ 1 — z. 



the n-particle transition amplitude can be written as 

A^-\K, h„ {h^}) E E b, c)S^f^'^ {z)At'\K, {hi}) . (B.7) 

de he 

Although the structure of this equation remains valid at all orders in perturbation 
theory, we will restrict in the following to the leading-order splitting. Therefore, 
C{de,b,c) is the color factor of the ggg vertex, Sgg^''^^ is the leading-order splitting 
function for fixed helicities (which can be obtained from the entries of table |I| as 
explained in the caption) and wA^""^"* is the transition amplitude for the process 1+2 
X -\- g, when the emitted gluon has color d^. The dependence upon the color indices 
{di} of the remaining n — 2 particles is not explicitly indicated in A^"'~^\ In the 
following, we will systematically suppress the indication of the dependence upon all 
color labels with the exception of de, which is the only one relevant for our purposes. 
The quantity {hi} is the set of the helicities of the n — 2 particles which are not 
involved in the splitting process. Notice that we only consider the range z < 1, where 



the splitting functions with fixed helicities are meaningful. Squaring eq. (|B.7|) and 
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summing over the color of the emitted gluon we get 



b,b' c,c' de,d'^ he,h'^ 



X At'\he,{hi}){A^-'\K,{hi})y . 



(B.8) 



Since C(a, 6, c) = fabc, and using the normahzation conventions of Mangano and 
Parke iTSl, we have 



E E Sbb'Scc'C{d,, b, c)C*{d'^, b', c') = 2CJd. 

b,b' c,c' 



(B.9) 



Therefore, eq. (|B.8|) becomes 



de \ he 



T.s^'t''''i^)^t'\K,{hi})] ■ (B.IO) 



Writing exphcitly the sums over he and h'^ and carrying out the multiphcation we 
obtain 



h„ {hi}) ' ^ M^^\h,, hj, {hi}) ' + 7^(/^„ h„ {hi}) 



where 



Af(-\h,,h„{hi})^ = 2glC, S^^'^'^{z)rY.\^t'\+dhi}) 



+ 2glC. S-^^^^{z)\Y.\^t'\-Ahi}) 



nih,,h„{hi}) = AglC,Rel^S^^'^^^iz){S;^^''^iz)y 
X T.^V\+dhi}){At'\-,{hi}) 



By construction, given eq. ( |B.7D , the quantity 



A^"-'\±,{hi}) =EEKr ^(±^{M) 

{di} de 



(B.ll) 



(B.12) 



(B.13) 



(B.14) 
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is the {n — l)-particle amplitude squared and summed over colors, for a given helicity 
configuration {±,{hi}). Therefore, taking into account 

1 1 



99 



(B.15) 



where P^g is the unpolarized Altarelli-Parisi kernel for z < 1, from eq. ( |B.12| ) we get 



hi,hj,{ht} {di} 



where 



h,{hi} 



(B.16) 



[B.17) 



is the full {n — l)-particle transition amplitude squared and summed over helicities 
and colors. From it, we define the Born amplitude squared by multiplying by the flux 
factor and averaging over the spin and color of the incoming particles 



_A/f{«-i) 



1 



fB.18) 



2ki ■ k2 uj{ai) uj{a2) 

When considering three-jet production at the next-to-leading order, this quantity 
coincides with A4^^'^\ defined in eq. ( |2.10| ). Eq. ( |B.13[ ) can be simplified using the 
explicit form of the splitting functions with fixed helicities. We have 

z{l-z) 



ol2 



(B.19) 



Therefore we get 



J2'^{hi,hj,{hi}) = 327casC^z{l-z) 

hi,hj,{hi} {di} 



X 



U*.yj {hi}d„{di} 



-,{hi}) 



We can now write the full n-particle amplitude squared 



(B.20) 



(B.21) 



2ki ■ k2 uj{ai) uj{a2) 

For three-jet production at the next-to-leading order, Ai^"^^ coincides with Ai^^\ 
defined in eq. Using eqs. and j ^^Uj ) we get 

IGTia.c 



C^z{l-z)M 



(n-l) 



B.22 



-41- 



where 



1 



1 



X Re 



2ki ■ k2 uj{ai) uj{a2 
' <i]> 



y-Ji {hi}d,,{di} 



(B.23) 



Notice that, apart from the factor < ij > /[ij], TW^"^^-' would be equal to the Born 
amplitude squared, eq. ( |B.18|) , if the first entry of A^"^^^ were equal to the first entry 
of (^(""^))*. In eqs. ( |B.18| ), ( |B.21| )-( |B^ ) we suppressed the obvious functional de- 
pendence, but we emphasize that, while A^^""^) is independent from fcj and kj, A^("~^) 
has a residual dependence upon these momenta through the factor < ij > /[ij], which 
explicitly appears in eq. ( [B.23| ). 



We now turn to case b), that is the splitting g — * qq. We use again eq. ( |B.6|) to 
define the z parameter, and assign color indices k and k to the quark i and to the 
antiquark j respectively. Eq. ( |B.7| ) still holds, with the formal substitution 



C(4,6,c) C{de,k,k) = (A'^^) 



kk 



(B.24) 



Obviously, we have also to substitute everywhere P^g with P^. Instead of eq. ( p.9| ) 
we have 

E E Skk'S-k-k,C{de, k, k)C*{d'^, k', k') = 2T,5,,4 . (B.25) 
fc.fc' k,k' 

The derivation goes unchanged from this point on. We have only to take into account 
that 

2 1 1 ^ 



E 1^^^'^' 

hi,hj 



E 1^- 

hi,h 



and 



Qg 



z(l-z) 



'hihj 



(B.26) 
(B.27) 



Eq. ( p.22| ) gets therefore modified as follows 



' kj 



+ l^T,z(l-.)A<("-) 



(B.28) 



Finally, we deal with the splitting q —>■ qg- With obvious modifications due to the 
fact that the exchanged particle in now a quark instead of a gluon, we can follow the 
derivation outlined before. Nevertheless, we immediately understand that in this case 
there is no TZ term. In fact, we know from eq. ( |B.10[ ) that such a term is present if and 
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only if for a given external helicity configuration the exchanged virtual particle can 
have both helicities. In the present case, this is forbidden due to helicity conservation 
along the quark line. More formally, for the q qg splitting we have always 

S^'-'' (^) (^)) * = St'-'' (^) {S;'-"' (^)) * = (B.29) 

for every possible choice of hi and hj. 

We collect now the results obtained so far, also explicitly indicating the depen- 
dence upon the momenta of the particles, referring to eq. (p3.6| ) for the definition of 
the kinematics: 



{ki,k2; ..,ki,..,kj,..) 



^"""^ ^a<5K,a,)W-M^"-'^ {kuk2;..,k,,..) 



U U aiS{ai,aj)^ 



+ ^ Qa.Sia.,a,Y{z)M^--'^ {k,, k,' .., fc,, .., fc,, ..) , (B.30) 

fv'i ' rCj 



where 



Qgg.{z) = -AC^zil-z), (B.31) 

Qgg.iz) = AT^zil-z), (B.32) 

Qgg*{z) = 0, (B.33) 

Qg,*{z) = 0. (B.34) 

In these equations, the * symbol over the flavour of the particle that eventually splits 
reminds that this particle is off-shell. In principle, this notation should be extended 
also to the Altarelli-Parisi splitting kernels, but at the leading order Pat* = Pa*b, and 
therefore there is no need to keep track of the off-shell particle. By construction, the 
A term of section 4 is equal to Q A^*-""^^. 

We have now to deal with the collinear emission from an incoming particle. We 
have two options: we can perform an explicit calculation, as we have done for the 
final state emission; or otherwise we can exploit the crossing symmetry property of 
the n-particle amplitude We will pursue this second option, and to be definite 



we will restrict to emission from particle 1. To proceed further, we have to think 
in terms of an (unphysical) amplitude, for which also the particle 1 is outgoing; the 
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limit 1 II J will therefore be given by eq. (|B.30| ) (with obvious modifications in the 
notation). The kinematics will be given by eq. 



ki = zkp , kj = {1 — z)kp , (B.35) 

where ki = —ki since in the unphysical amplitude the particle 1 has negative energy. 
In doing the crossing, apart from the substitution ki ^ ki, we will have to refer the 
kinematics to ki, that is 

kj = {1 — z)ki , kp = zk\ (B.36) 

or 

z — ^ l/z. (B.37) 

Since kp is the four-momentum of an incoming particle in the {n — l)-particle ampli- 
tude, the second equation in ( |B.36[ ) implies that the flux factor of the {n — l)-particle 
amplitude will be I/2; times the flux factor of the n-particle amplitude. Therefore, to 
have the correct normalization, we have to multiply by z the (ra— l)-particle amplitude. 
Furthermore, since the flavour of 1 and P can be different, a factor uj{ap) / uj{ai) will 
take into account the different color and spin normalization factors of the n-particle 
and {n — l)-particle amplitudes. Eq. ( p.30|) becomes therefore 



(fci,A;2; ..,A;„..)^ 

_ i^^^^^^^^^ii^^P< (-\M^^-^Uzk. k.- 1 

^ ^ fci-fc, u:{a,) '^-^^sin,,a,){^)M {zkuk,,..} 

ki ■ kj uj{ai) ^ ^' \zj 



(B.38) 



where a{g) = and cr{q) = 1 (therefore, (— takes into account the crossing of a 
fermionic line), and the overall minus sign is due to the ki — > ki substitution. We can 
now exploit the crossing symmetry property of the Altarelli-Parisi kernels (for z < 1) 



uj{a) \zj 

For consistency with this equation, we are thus led to define (notice the * symbol) 



Qb^aiz) = -(-)K'^i+'^^''^J ^ - . (B.40) 

uj{a) 
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Eq. ( [B.38|) becomes therefore 

+ Qs(a,,-a,ra, iz)M^--'^ (zk,, k^. .., fc,, ..) , (B.41) 



where, from the definition in eq. ( |B.40|) , and using eqs. (|B.31|) -( p.34]) 



Qy.y{z) = -4C^l^, (B.42) 
Q,.g{z) = 0, (B.43) 

Qg.,{z) = -4C,^, (B.44) 
Q,*,{z) = 0. (B.45) 

In spite of the fact that the A term has no deep physical meaning, still the crossing 
symmetry property has to hold true, and the relation between Qab* and Qb*a should 
not be regarded as purely incidental, that is, due to the definition ( |B.40|) . In fact, 
directly performing the calculation for the initial state collinear splitting, and defining 
Qh*a through eq. ( |B.41[ ), we obtain again eqs. ( |B.42| )-( p.45| ). 



We finally list the results for Ai^^') for processes in which only massless partons 
are involved. This kind of processes are relevant for three-jet production. The corre- 
sponding Born amplitudes squared, are well known from the literature (see for 
example ref. ||30||). 

We start with the purely gluonic processes, denoting the momenta of the glu- 
ons which do not split by k,l,m,n, and with P the momentum of the gluon which 
eventually splits. We get 



1 AglN^jN! - 1) 
2ki ■ k2 uj{ai)uj{a2) 

X Re |M ( lkl)Hmn)^ + (km)Hln)^ + (kn)UlmY 

[ m ^ 

X E (^{kl){lm){mn){nP){Pk)y^\. (B.46) 

a*{k,l,m,n) J 

Here a*{k, 1, 171,12) denotes all the permutations of the elements k,l,m,n which are 
inequivalent under reflection {{klmn} {nmlk}). For the processes with two quarks 
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and three gluons we denote by m and n the momenta of the gluons which do not 
participate in the sphtting, and by q and q the momenta of the quark- antiquark pair. 
We have 



2ki ■ k2 u;{ai)u;{a2) 



X Re 



X 



{Nl - If 



({qm){mn){nP){Pq 



'^^c aim,n,P) 

^ ..(i;,P) [{qm){mn){nP){Pq){qn)^{mP)^ + ^ 



a* (m,n,P) (^^) ("^^) (^-^) (-^^) ("^^) '^{Pri)^ 



(B.47) 



Finally, for processes with four quarks and one gluon, we have to distinguish between 
the case in which there are two quark- antiquark pairs of different flavour, and the 
case in which the flavour of the two pairs is the same. In both cases, we denote the 
momenta of the quarks by q, q, Q and Q. The momentum of the splitting gluon is 
again P. For quark-antiquark pairs of different flavour we get 

2h ■ k2 u{ai)u;ia2) \ [ij] {qq?{QQ? 



{qQ? , {qQ) 



2 



+ 



{QPYiqPy {qPY{QP)\ 
Nl-lf {qq? ^{qq){qQ) 2 {qq){qQ) 



K \{qP)HqPy {qP)^{qP){QP) {qP)^qP){QP) 

jQQY _ _2 {qQ){QQ) _ _jmmi_\ 

^ {QP)'{QP)' {QP)HqP){QP) {QP)^qP){QP)J 
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For quark- ant iquark pairs of equal flavour we have 



1 4gl{N^-l) 



2ki ■ k2 uj(ai)uj(a2) 



m {qq){qQ){qQ){QQ) 



X 



{qqy 



{qQY 



+ -7^ 



{Qqf 



+ -7^ 



{qPy{qPy {qP?{QP? {QP?{qP? {QP)HQP) 



{qq){qQ) 



+ -7^ 



miQQ) 



iV2 \{qP)^qP){QP) {QPy{qP){QP) 



{qq){qQ) 



{Qq){QQ)_ \ 

{qPymiQP) {QPy{qP){QP)J 



[B.49) 
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